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We report the prediction of quasi-bound states (resonant states with very long lifetimes) that 
occur in the eigenvalue continuum of propagating states for a wide region of parameter space. 
These quasi-bound states are generated in a quantum wire with two channels and an adatom, when 
the energy bands of the two channels overlap. A would-be bound state that lays just below the 
upper energy band is slightly destabilized by the lower energy band and thereby becomes a resonant 
state with a very long lifetime (a second QBIC lays above the lower energy band). 



Since the bound state in continuum (BIG) was first 
proposed in 1929 by von Neumann and Wiener [li, vari- 
ous researchers have reported its existence [1, [3,13, S 0, 
0, H 0] ■ All studies agree that the phenomenon can only 
occur at discrete points of parameter space (i.e., the BIG 
is a zero- measure effect). 

We here report the existence of a gwasz-bound state in 
continuum that exists over a wide region of parameter 
space (finite measure). By quasi-bound state, we mean 
a resonant state with a very long lifetime, so long that it 
appears to be a bound state in space and hardly decays 
in time. The quasi-bound state emerges when the system 
has an impurity level coupled to two overlapping energy 
bands with divergent van Hove singularities at the band 
edges. If there were just one energy band, a bound state 
would appear just outside the edge (due to the singular- 
ity). This would-be bound state is slightly destabilized 
by the other energy band, forming a quasi-bound state. 
Since the quasi-bound state is virtually a bound state 
in continuum, it may be useful as a high-energy excited 
level of a carrier as a quantum device application. 

The Hamiltonian for our system is a tight-binding 
model on a ladder with an adatom (Fig.[lja)): 
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Here, cj. y and Cx,y are the creation and annihilation op- 
erators of a spinless fermion at the site (x, y) with integer 
X (—00 < a; < 00) and y = 1,2, whereas and d repre- 
sent a dot, or an adatom, attached to the (0, 1) site of the 
ladder. The first line of Eq. ([T]) gives the hopping matrix 
elements along the ladder, the second line the hopping 
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FIG. 1: (a) An adatom, or a dot, attached to a ladder, (b) 
After diagonahzing the ladder, the system is composed of a 
dot coupled to two independent channels. 



elements across the ladder and the third line gives the 
hopping elements to and from the one-particle level of 
the dot. 

The ladder has two eigenmodes Cx,± = icx,i±Cx,2)/ V^, 
which transform the Hamiltonian JT]) to 
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see Fig. [Ijb). The system has two conduction channels 
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with — TT < fc± < TT, each of which forms an energy band 
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Two divergent van Hove singularities occur in the den- 
sity of states function at the edges of both bands 14 1. 
The two bands overlap (with one of the singularities of 
one band embedded within the continuum of the other) 
whenever \t'-^\ < |th|- The dot level couples to both chan- 
nels as can be seen in Eq. Incidentally, we can also 
regard the Hamiltonian ^ as conducting electrons with 
spin cr on a tight-binding chain under a magnetic field 
proportional to t'-^. 

In terms of the new operators for the wave number k 
in Eq. ([3]) , the Hamiltonian ([2]) takes the form of the cou- 
pled Friedrichs-Fano (Newns- Anderson) model Hamilto- 
nian [l^, 11, 13, 13]) which has been thoroughly studied. 
Using the standard argument for this model, we obtain 
the dispersion equation for this coupled system: 
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This is equivalent to a twelfth order polynomial equation 
for z. The complex solutions of this equation correspond 
to the complex energy eigenvalues of the resonance states, 
with the decay rate given by the imaginary part. 

Here, we focus on the numerical solution of this equa- 
tion, using = 0.345th and g — 0.1th for our demonstra- 
tion. The twelve eigenvalues for Ed — 0.3th are listed in 
Table I, as an example. Each discrete eigenvalue can be 
distinguished by its position on the complex surface, 
the complex K- surface and the complex energy surface, 
where the three quantities are related by 
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The corresponding eigenfunction is given in the form 
f *(x,l;t)\ ^ 
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with appropriate constants A±. 

The complex energy surface is composed of four Rie- 
mann sheets. It is a physical requirement [l5| that a 
resonant state decays in time but diverges in space. The 
resonant state, therefore, has a negative imaginary part 
for its energy and a negative imaginary part for one or 
both of the wave numbers if + and . An eigenvalue 
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FIG. 2: (a) The wave function modulus \'il)(x,y)\ of the state 
Q2 around the origin, (b) The same but away from the origin 
on the logarithmic scale. The plots for y = 1 (the upper 
leg) and y = 2 (the lower leg) are almost indiscernible. The 
parameters are set to t{-^ = 0.345th, g ~ 0.1th and Ed = 0.3th. 



in the upper plane and in the upper plane is de- 
fined to lay on Riemann sheet I and is a complete bound 
state; the states PI and P2 in Table [J are such states. 
An eigenvalue in the lower plane and in the lower 
K- plane is defined to lay on Riemann sheet IV and is a 
complete resonant state; the state SI in Table U is such a 
state. (A state with a positive imaginary part of the en- 
ergy is a so-called anti-resonant state, which is the time 
reversal of a resonant state and hence grows in time; S2 
is the time reversal state of SI.) 

An eigenvalue in the lower plane but in the upper 
K- plane is defined to lay on Riemann sheet II, while 
one in the upper plane but in the lower plane is 
defined to lay on Riemann sheet HI. The resonant states 
on these two sheets (states Q2, Q4, and R2 in Table |l| 
are the main focus of the present paper. 

The resonant states in sheets II and HI are resonant 
states because they diverge in space due to the negative 
imaginary part of only one of the wave numbers; see Eq. 
([7]). Some of them, however, diverge in space very slowly 
and decay in time very slowly (the latter is demonstrated 
in Table [l|. The state Q2, for example, appears to be a 
localized state around the x-axis origin (Fig. [H^a)) and 
exponentially diverges only far away (Fig. [2Ib)). Indeed, 
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TABLE I: The twelve discrete eigenvalues for 4 = 0.345th, g = 0.1th and Ed = 0.3th. 



state 


E/th 


K+ 


K- 


Riemann Sheet 


Pl(0.3) 


1.34501152 






3.14159265 


+i 1.11593256 


3.14159265 


+i 0.00480148 
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P2(0.3) 


-1.34500463 








+i 0.00304629 




+i 1.11592751 


I 


Ql(0.3) 


1.34501136 






3.14159265 


-i 1.11593245 


3.14159265 


+i 0.00476787 


II 


Q2(0.3) 


-0.65501370 


~i 1.5093 xlO" 
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1.25558888 


-i 1.5875 xlO"^ 


-0.00002882 


+i 0.00523534 
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-0.65501370 
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0.00002882 


+i 0.00523534 


II 


Q4(0.3) 


0.29998854 


-i 0.00153774 




2.27180290 


-i 0.00201224 


-1.52576970 


+i 0.00153930 
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Rl(0.3) 


-1.34500459 








+i 0.00303273 




-i 1.11592748 


III 


R2(0.3) 


0.65509906 


~i 2.9331 xlO" 
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III 


R3(0.3) 


0.65509906 
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-i 3.0852 xlO"'* 


III 


Sl(0.3) 


0.29991927 


-i 0.01154476 




2.27161773 


-i 0.01510419 


1.52570333 


-i 0.01155625 


IV 


S2(0.3) 


0.29991927 +i 0.01154476 
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-i 0.01510419 


-1.52570333 
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solving the dispersion relation ([5]) by perturbation ex- 
pansion in we can show that the imaginary part of the 
energy of state Q2 is proportional to (extremely small 
for g I) due to the interaction between the divergent 
van Hove singularity in the upper energy band and the 
continuum of the lower band. Since the real part of the 
energy lays within the lower energy band, we refer to this 
state as a quasi-bound state in continuum (QBIC). 

In particular, the quasi-bound state Q2 has the real 
part of the energy just below the lower edge of the upper 
energy band e_, which occurs at — + = —0.655th in 
the present case. If we had only the upper channel, this 
state would be a bound state (due to the van Hove sin- 
gularity) In order to elucidate how the QBIC effect 
appears, we compare the above result to that obtained 
for the separate one-channel systems with the following 
new Hamiltonians with a structure similar to Eq. ([2|): 

a; — — oo 

+ ^ (d''cQ + 4d) + EddU T 4- (8) 

Here we set the coupling between the site x = and the 
adatom to g/V^ for quantitative comparison with the 
Hamiltonian ([2]). We have also added the energy offset 
=Fih so that the energy band of each Hamiltonian Tl± 
may coincide with the energy bands e± of the ladder 
system. The Hamiltonian Ti- thus mimics the upper 
energy band £_ of the — channel of the ladder system, 
while the Hamiltonian 7i+ mimics the lower band e+. 

The bound state below the lower band edge of the 
Hamiltonian for t[^ ~ 0.345th, g = 0.1th and = 
0.3th has a pure real energy E/th ~ —0.65501371 and 
a pure imaginary wave number K = iO. 00523550. This 
bound state of the one-channel system (Fig. (Ha)) indeed 
closely resembles the eigenvalue and the wave number 
of state Q2 of the two-channel system (see Table HI. 



When the second channel is attached to the one-channel 
system (Fig. |31Ib)), the energy is inside the conduction 
band of the -I- channel so that a portion of the bound 
particles outside the edge of the — channel leak into the 
-I- channel and escape to infinity. Notice KeK- < 0, 
while KeK+ > 0; the particles on the channel — are suck 
into the origin and spring out of the origin onto the -I- 
channel. This leak makes the state Q2 a resonant state; 
it can be generally shown [l5| that the imaginary part of 
the resonant eigenvalue arises due to the momentum flux 
escaping from the scattering potential. 

The actual eigenfunction of the state Q2 is exemplified 
in Fig. [3jc). One can show that the amplitude at the 
dot tp{d) is order g smaller (g = 0.1th in this case) than 
the amplitude of ■!/'(0, — ) and the amplitude of 'ip{0, +) is 
order g smaller than tp{d). Hence the amount of the leak 
is small (^ g^). The same analysis applies to the state 
R2, which is also a quasi-bound state. Figure [4] shows 
that the quasi-bound states survive over a wide range of 
the parameter E^. 

To summarize, we have found quasi-bound states in 
the continuum over a wide range of system parameters, 
which is in striking contrast to the bound state in the 
continuum. The quasi-bound states emerge in the system 
with two overlapping eigenvalue continua when at least 
one of the bands contains a divergent Van Hove singular- 
ity. For instance, a would-be bound state just below the 
upper continuum turns out to be a quasi-bound state in 
the lower continuum. The simple picture of the mecha- 
nism suggests that the quasi-bound state may be found 
in various systems including a microwave tube with an 
ion [3] and in the scattering of ions and nuclides. 
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FIG. 3: (a) A schematic view of the strongly bound state 
(due to the van Hove singularity) of a one-channel system 
with the eigenvalue just below the lower band edge, (b) Some 
of the bound particles leak into the attached channel, (c) 
The eigenfunction of the state Q2 for t'y^ = 0.345ih, g = 0.1th 
and Ed ~ 0.3th- The amplitude modulus of the — channel, 
|'(/'(a;, — )|, that of the -I- channel, \^{x,+)\, and that of the 
dot, |V'(d)|! are indicated. 
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